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The properties of nearby galaxy clusters limit the range of cosmological parameters 
consistent with our universe. We describe the limits which arise from studies of the 
intracluster medium (ICM) mass fraction ficM 3,nd consideration of the possible 
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sources of systematic error: £7jv/ < at 95% confidence. We emphasize 

that independent of Type la supernovae (SNe la) observations, this cluster study, 
taken together with published cosmic microwave background (CMB) anisotropy 
studies, indicates a non-zero quintessence or dark energy component Oq > 0. 
We then discuss future galaxy cluster surveys which will probe the abundance of 
galaxy clusters to intermediate and high redshift. We investigate the sensitivity of 
these surveys to the cosmological density parameter Q,m and the equation of state 
parameter w of any quintessence component. In particular, we show that cluster 
survey constraints from a proposed large solid angle X-ray survey are comparable in 
precision and complementary in nature to constraints expected from future CMB 
anisotropy and SNe la studies. 



1 Overview 

Galaxy clusters contain a wealth of accessible cosmological information. Their 
masses range from 10^^-10^^ Mq, making clusters the most massive collapsed 
or virialized objects in the universe. Although numerous careful and varied 
studies of nearljp/.rliisters indicate that they are still accreting mass at the 
present epochtZHl2Hl3'E3, it is evident that clusters exhihitpStriking regularity in 
scaling relations between mass, size and temperaturec3'E3E3. In fact the scatter 
of clusters about the X-ray size-temperature relation is 159L-comparable to the 
scatter of elliptical galaxies around the fundamental planecJ. 

Both these characteristics (accretion at the present epoch and tight scaling 
relatioiis)|_are also exhibited by clusters formed within hydrodynamical simu- 
lations oEjtl, providing some confidence that the process of cluster formation 
is dominated by gravity and gas dynamics, and therefore simple enough to 
to be effectively modeled in numerical simulations E3. Ongoing observations 
with new X-ray observatories and proposed radio observatories will ena 
to further improve our understanding of cluster formation and evolution 
Because galaxy clusters are sufficiently regular that their masses can be esti- 
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mated from observables such as the ICM temperature T^, yields from cluster 
surveys of the high redshift universe are more readily interpreted and can. in 
principle, be used to constrain cosmological parameters. 



2 Constraints on VLm from ICM Mass Fractions 

Galaxy clusters can be used to study the mix of baryonic and dark matter on 
scales of roughly 10 Mpc. Because there are no candidate mechanisms available 
to segregate baryons and dark matter on these scales E3, it is often argued that 
the the baryon fraction within clusters fd should reflect the universal baryon 
fraction fs = VLb/^m, where ^Ib {^m) is the cosmological density parameter 
of baryons (all clustered matter). Therefore, a measure of the cluster baryon 
fraction can be combined with primordial nucleosynthesis constraints on 
the baryon to photon ratio and measurements of the CMB temperalure to 
yield an estimate of the cosmological density parameter VIm = ^b/ fcV^- 

There are at least three reservoirs of baryons in galaxy clusters: (1) the 
ICM, (2) the galaxies, and (3) dark baryons. The X-ray bremsstrahlung and 
recombination radiation from the ICM provides a precision tool for estimating 
the baryonic mass in the ICM reservoir. The optical light emitted by stars 
within galaxies provides a more blunt estimate of the baryonic mass in galaxies 
(requires accurate estimates of the baryonic mass to light ratio in typical cluster 
galaxies), and to date there is no clean way of separating a possible dark 
baryonic component from the dominant dark matter reservoir within clusters. 
Detailed studies of individual clusters teruLto indicate that the ICM mass is 
several times larger than the galaxy massll3.Below we describe a study which 
uses observations of the ICM baryon reservoir to place an upper limit on ^m- 



2.1 Study of the ICM in an X-ray Flux Limited Sample of 45 Clusters 

The ICM density profile p{r) can be extracted from an image of the cluster X- 
ray emission, given an estimate of the mean ICM temperature , an emission 
model and some assumption about the cluster geometry. The central X-ray 
surface brightness can be expressed 

= ^77tVi4 r ^^HT.), (1) 

where A is an emission coefficient, z is redshift, vrip is the proton rest mass, 
Ui = pI ^ii^ and rii is the number density of species i. We use ROSAT PSPC 
obseijiiations of 45 clusters from an X-ray flux limited sample with available 
dataB Within our chosen PSPC band (0.5:2 keV), the emission coefficient 
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A is extremely insensitive to ICM temperature CJ, and so one can accm'ately 
infer p{r) without direct knowledge of the ICM temperature structure. One 
does need to assume a geometry; we assume the clusters are spherical. This 
introduces errors in the density inversion, and we characterize these errors by 
testing our analysis on mock observations of simulated galaxy clusters. Using 
mock observations of 48 hydrodynamical simulations, we calculate that the 
cluster ICM mass within a radius 7f§flOj roughly half the virial radius, can be 
estimated with an accuracy of 10% Ej. 

Figure 1: Measured ICM mass fractions 
fiCM versus mean ICM temperature Tx 
for an X-ray flux limited sample of clus- 
ters. The mean ficM for the clusters with 
ksTx > 5 keV (vertical line) is 0.212h~^^'^ 
(horizontal line). This measurement pro- 
vides a lower limit on the fraction of matter 
within cluster virial regions which is bary- 
onic. Together with current best estimates 
of the scale of systematic errors on this up- 
per limit, this measure of (fichl) provides a 
95% confidence upper limit on the cosmolog- 
ical density parameter of clustered matter: 




Qui < OAAh', 
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To estimate the ICM mass fraction ficM we not only need the ICM den- 
sity profile p{r), but we also need the cluster binding mass. We estimate M500, 
the binding mass within rsoo, by assuming the ICM is in hydrostatic equilib- 
rium and is isothermal. Departures from equilibrium and isothermality will 
introduce errors. A common temperature profile in clusters could potentially 
lead to sys tem atic errors in our binding mass estimates. These are considered 
below in J 2.2 . 

Fig. |l| contains a plot of our ficM measurements versus emission weighted 
mean ICM temperature in 45 clusters. There are two important charac- 
teristics of the distribution of ficM- First, there is a weak, but statistically 
significant, tendency for low mass (low ksTx) clusters to have lower ficM- The 
physics responsible for depleting the ICM in low mass clusters is thought to be 
preheating of the intergalactic medium by star formation within_galaxies be- 
fore the gas collapsed into the forming potential wells of clustersBEj. Generally 
speaking, preheating of the gas prior to cluster formation or energy injection 
after cluster formation will have a larger effect on low mass clusters than on 
high. Thus, in using ficM to constrain the density of clustered matter, we 
restrict ourselves to the highest mass, hottest systems: ksTx > 5 keV. 

Second, splitting the sample at 5 keV (vertical line), we find the mean 
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fiCM = (0.212 ± 0.006)ft.5(f/^ (the horizontal line) . J:|lii!i,numbcr IS m reason- 
ably good agreement with other estimates of //cA^liy'cl. Constraints on the 
primordial deuterium abundance from high redshift absorption ssffitems pro- 
vide an estimate of the baryon density parameter: Qb = 0.076/i5oQ. Together 
with our lower limit on the cluster baryon fraction /cj, this number implies an 
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upper limit f^M < (0.36 ± 0.01)/i5o , where the uncertainty is only statistical. 
We require a factor of three systematic error to reconcile this number with 
Hm = 1 cosmological models. 

2.2 Discussion of Systematics 

It is particularly interesting to consider systematics which could potentially 
make f2j\/ = 1 models more consistent with our data. There are several possible 
systematics: (1) any mechanism which enhances the cluster baryon fraction 
relative to the universal baryon fraction, (2) overestimating the ICM mass, 
and (3) underestimating the cluster binding mass M^qq. 

First, there are no known mechanisms for enhancing the baryon fraction by 
a factor of three on scales of 10 MpcE3. Second, plausible effects which would 
cause us to overestimate the ICM mass are (a) the effective area of the ROSAT 
PSPC is known to only 15%, corresponding to a 7.5% systematic uncertainty 
in our measured ICM masses, and (b) clumping or multiphase structure in the 
ICM would enhance the X-ray emission relative to single phase gas. In fact, 
our hydro simulations indicate that th£_X-ray emission is enhanced by roughly 
20% on average by clumping in the gasl23; this leads us to overestimate the ICM 
mass by 10%, on average (we have already corrected for this). Interestingly, 
preliminary results from an ongoing analysis of an independent set of hydro 
simulations indicates that clumping from infalli|ag substructure is present at 
about the same level as seen in our simulations In addition, this clumping 
can potentially be addressed by comparison of ICM mass fractions derived 
from X-ray data and those derived from Sunyaev-Zel'dovich Effect (SZE) ob- 
servations. Because the SZE is sensitive to a line integral of the pTe, where 
is the electron temperature, clumping would likely have a much stronger effect 
on the X-ray measures than on the SZE. A comparison of ficM derived from 
an analysis of SZE observations of 18 clusters provides no indication for sys- 
teraaiically different results within an accepted range of the Hubble parameter 
Horn. 

Third, there are many possible effects which could systematically affect 
our binding mass estimates M500; these include bulk flow or turbulence, sup- 
port from magnetic fields, and ICM temperature profiles. Indeed, analyses of 
some individual clusters provides binding masses derived from galaxy dynamics 
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which are significantly higher than those derived from hydrostatic equihbrium. 
However, a systematic study by the CNOC cohaboration of 14 intermediate 
redshift clusters finds tha^ratio of galaxy dynamical to isothermal hydrostatic 
masses to be 1.04±0.07l3. Further studies using gravitational lensing and the 
spatially resolved ICM temperatures available from Chandra should provide 
much needed additional information. 

Taking the 7.5% systematic from the uncertainty in the PSPC effective area 
and the 7% systematic uncertainty on the binding mass estimates, we estimate 
a total systematic uncertainty of 10%. This together with the observations 
outlined in the previous section leads to a 95% confidence upper limit of JIa/ < 

2.3 Two Independent Observational Arguments Jar 51q > 

Studies of high redshift SNe la prefer cosmological models with Via > 
Both cluster baryon fraction arguments and mass to light studies u favor low 
U,M models. The mass to light ratio studies are more difficult to interpret, 
because the stellar populations of galaxies inside clusters differ significantly 
from those outside clusters. Nevertheless, these two approaches, subject to 
different sjzatematics, indicate Q,m « 1- Together with constraints on CMB 
anisotropjnil, these clusters lead to the conclusion fig > 0, independent of the 
SNe la studies. 



3 Cluster Surveys and Cosmology 

The relatively simple evolution (compared to galaxies) and regularity of galaxy 
clusters make them candidate tracer particles to_use in measuring the volume- 
redshift relation. This classical cosmological testLJ has been applied to galaxies 
with limited success, due at least in part to the complex relation between galaxy 
brightpfiss and mass and the poorly understood evolution of the galaxy abun- 
danceE3(but see Newman & Davis article for discussion of new approach being 
considered in the DEEP survey). Clusters are more amenable to these studies, 
because our theoretical understanding of their structure and evolution is more 
complete. We don't expect the abundance of clusters to remain constant with 
redshift, but we can calculate its evolution, enabling the volume-redshift re- 
lation test. Moreover, the cosmological sensitivity of the abundance evolution 
itself provides additional leverage. 

Galaxy cluster surveys of the nearby universe are an old endeavorEI; how- 
ever, new technology and techniques are now making it possible to carry out 
extensive cluster studies of the intermediate and high redshift universe. Stud- 
ies of particular interest include: (1) a 400 deg^ serendipitous XMM cluster 
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survey extending to high rcdshift Ell, (2) the 10'* deg^ survey of the nearby 
iand intermediate redshift universe with the Sloan Digital Sky Survey (SI 
□, (3) a 12 deg^ interferometric SZE survey of the high redshift universeC 
and (4) a 10^ deg^, deep X-ray survey of the nearby and intermediate redshift 
universe. 

Figure 2: We plot the expected redshift 
distribution of a proposed SZE survey for 
three different cosmological models. The in- 
set (upper right) shows the mass of clus- 
ters detectable at 5(7 significance as a func- 
tion of redshift. Note the mild redshift sen- 
sitivity of this mass threshold; this is the 
unique characteristic of an SZE survey. This 
particular proposed interferometric SZE sur- 
vey wiW cover 12 deg^ in one year, and it 
should yield a sample of roughly 400 clusters 
if the currently favored cosmological model 
(f2j\/ = 0.3, Ca = 0.7, erg = 1) is correct. 
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Fig 1^ contains expectations for the redshift distribution of the interfero- 
metric SZE survey proposed by J. Carlstrom and collaborators. Because the 
SZE is a distortion of the CMB spectrum caused by inverse Compton scat- 
tering of CMB photons with hot electrons in the ICM, the SZE doesn't suffer 
from the cosmological dimming that any light source experiences. This fact 
makes the SZE ideally suited to studies of massive structures in the high red- 
shift universe. The inset of Fig ^ contains the mass of a_cluster which would 
be detected with 5a significance as a function of redshifca. Interestingly, this 
limiting mass is relatively constant with redshift, and it corresponds to a very 
low mass galaxy cluster; thus, the proposed Carlstrom survey will enable us 
to probe the universe for clusters to the very moment of their emergence. As 
described below, this fundamental observation will provide a powerful test of 
structure formation models and allow precision measurements of several cos- 
mological parameters. 



3. 1 Constraining the Equation of State Parameter w 

Because there are now two independent observational arguments for a non-zero 
quintessence or dark energy component Qq, further work is required not only 
to test this conclusion, but also to make measurements of the equation of state 
parameter w = p/p of this component. In principle the proposed SZE survey 
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Figure 3: We show the sensitivity of an SZE survey to changing w (left) and changing Qm 
(right). These differences are separated (top) into the effects on the volume surveyed and 
the abundance evolution. Both the volume (light line) and abundance (heavy line) for each 
model are shown relative the the volume and abundance of the fiducial model {w = —I and 
= 0.3). Note the clear differences between changes in w and changes in Slj\f , indicating 
that in principle w and Qm can be determined simultaneously. 

can do just thi0, because the equation of state of the dark energy detemines 
how its energy density evolves pq oc where R is the scale factor. 

This evolution affects the expansion history of the universe, which is coupled 
to the volume-redshift relation and the growth rate of density perturbation£3. 

Fig m contains a plot of the ft m and w dependence of the yields of the 
SZE survey in the case that the mass limit is taken to be a constant Mum — 
2 X IQ^^h^Q Mq. Note (lower left) that increasing w from -1 (the A case) 
to -0.2 ( with fixed Q^m = 0.3) decreases the number of cluster expected at 
intermediate redshifts, but increases the number expected at high redshift. 
This is explained in the upper left panel; at low redshift, the larger surveyed 
volume of the w = —1 model is the dominant factor, whereas at higher redshift 
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the evolution of the abundance plays a larger — and eventually dominant — 
role. The right panel contains expected yields as a function of Qm (with fixed 
w = —1). Lower flM yields are larger, and they grow ever larger fractionally 
with redshift. This is essentially the result of the evolution of abundances, 
because the volumes probed are very similar in these models (upper right). 
The slower growth of density perturbations in low ^Im models means that 
the abundance changes more slowly, and we expect to see clusters to higher 
redshift. For more realistic tests which include the cosmological dependence 
of the limiting mass Mumiz), please see Haiman, Mohr & Holder (2000). 

3.2 Complementary, High Precision Cosmological Constraints 

Finally, we turn to consider the cosmological constraints possible from the large 
solid angle, deep X-ray survey proposed by G. Ricker, D. Lamb and collabora- 
tors. Although this survey extends only to a redshift z ~ 0.7, it will detect 10** 
clusters and provide temperature measurements for ~2,000 of these. These 
large numbers make for strikingly precise cosmological constraints, and having 
the and luminosity measurements for these clusters makes it possible to-test 
for unusual evolution models while constraining cosmological parameters Eil. 

Figure 4: Here are estimates of the power 
contained in a large solid angle X-ray clus- 
ter survey to measure ^.nd w simultane- 
ously. The dots mark the Icr, 2cr and 3(t con- 
fidence regions with which models can be dif- 
ferentiated from a fiducial model (Qm = 0-3, 
ui = — 1, (T8 = l). We show only the slice of 
parameter space with h = 0.66; all models 
have Qm + Qq = 1. For comparison, the 
dashed lines indicate the constraint from a 
1% measurement of the location of the first 
doppler peak in the CMB anisotropy spec- 
trum, and the dotted lines indicate the con- 
straints from a 1% measurement of the dis- 
tance to redshift 2 = 1. In combination with 
cither the CMB or SNe la studies, cluster 
surveys enhance the constraints on w. 
w 

Fig ^ provides an estimate of the extent to which the clusters from this 
proposed survey will allow one to simultaneously constrain w and Qm- In 
these models we take flM + — 1, and we require the local abundance of 
cluster above some mass threshold to be the same in each model. For this 
figure we assume some fiducial model: 57m = 0.3, Qq — 0.7, ag, — w = — 1 
and h = 0.66. We then use the total number of expected clusters and their 
redshift distribution to quantify the differences between any given model and 
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the fiducial model (redshifts for these clusters will be extracted from the SDSS 
imaging and spectroscopic survey). Note the tight, simultaneous constraints 
on w and ^Im- For comparison, we show (dashed) the constraints correspond- 
ing to a 1% measurement of the location of the first Doppler peak in the CMB 
anisotropy spectrum, and (dotted) a 1% measurement of the distance to red- 
shift z = 1. Note that the cluster survey does better than a 1% measurements 
to z — 1. In addition, the parameter degeneracies from our cluster study and 
these future CMB and SNe la constraints are roughly orthogonal, making the 
cluster survey complementary to either of the other studies. This orthogonality 
stems from the cluster survey sensitivity to abundance evolution. 

3.3 Systematic Effects and Non-standard Evolution 

Studies of the evolution of cluster structure are an important component of 
the effort to use cluster surveys to constrain cosmological parameters. The 
structure of cluster virial regions can affect the relation between the cluster 
virial mass and observables like X-ray luminosity, emission weighted ICM tem- 
perature and SZE flux. Very specific evolution models follow from theoretical 
work, and the observations required to test these models are now more readily 
available. Nevertheless, still uncertain physics, such as the effects of heating of 
the intergalactic medium before cluster formation, provides a potential source 
of systematic error in interpreting cluster surveys; we and others are currently 
studying these effects using hydrodynamical simulations. The enormous suc- 
cesses in modeling galaxy clusters to date and the vast array of data soon to 
be available on intermediate and high redshift clusters lead me to adopt an 
optimistic view; I suspect that through further enhancing the hydrodynami- 
cal cluster simulations, we can tackle the finer points of cluster evolution and 
cosmology simultaneously. 
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